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Abstract 

A systematic study of charged hadron multiplicities (n c h) at various collision energies is very much 
important in understanding the basic production mechanism of the hadrons in nucleus-nucleus 
collision experiments. Furthermore, the variations of n c h in nucleus-nucleus collisions with respect 
to the colliding energy and mass number can provide a potential probe for the formation of 
quark gluon plasma (QGP) in the laboratory. In this paper, we propose a phenomenological 



model based on the constituent quark-quark interactions to calculate the average multiplicity 



(n c /j) and pseudorapidity density at mid-rapidity ({dn c h/ drfj^o) of charged hadrons at various 
center-of-mass energies (a/sjvjv) for nucleus-nucleus (A — A) collisions. We first propose a new 
parametrization for and [dn c hj drjf? =G in p — p interactions based on some initial inputs 
which fit the experimental data very well. We further extend this parametrization by using 
simple phenomenological assumptions regarding mean number of participating quarks and mean 
number of collisions to obtain the n c h and (dn c h/ 'drj) rj=o for A — A collisions and show their 
dependencies on the mass number of colliding nuclei as well as on y / s/v r /v. We also compare the 
results obtained from our model with the results obtained from the modified Glauber model in 
order to demonstrate the difference between the two formalisms. Finally, we compare the charged 
hadron multiplicity and pseudorapidity density at mid-rapidity for A — A collisions obtained from 
our model with the available experimetal data from various heavy-ion collision experiments and 
give our predictions for A — A collisions at the Large Hadron Collider (LHC) and at Compressed 
Baryonic Matter (CBM) experiments. 
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I. INTRODUCTION 



The ultimate goal of ultra-relativistic heavy-ion collision experiments is to test the predic- 
tions of quantum chromodynamics (QCD), which is unanimously believed to be the theory of 
strong interactions [1]. The ultra-relativistic nuclear collisions reveal the nature of hadronic 
interactions at very short time and/or distance and throw light on the role played by the in- 
ternal structure of hadrons in the multiparticle production. Intensive theoretical efforts have 
also continued for more than three decades to understand the particle production mechanism 
in the nucleus-nucleus interactions. Hadronic multiplicities and their correlations can reveal 
information on the nature, composition, and size of the fireball from which they are origi- 
nating [2]. Of particular interest is the limit up to which the chemical/thermal equilibration 
is possible in the fireball. The appearance of quark gluon plasma (QGP) which involves a 
partonic medium at local thermodynamic equilibrium, and its subsequent hadronization af- 
ter the phase transition, should in general drive hadrons also towards chemical equilibrium. 
The theoretical description of particle production depends essentially on two ingredients : 
many-body interactions and the understanding of the production processes. 

The search for some regularities and systematics in the multiplicity distributions in 
hadron-hadron, hadron-nucleus and nucleus-nucleus collisions is a fascinating topic because 
it hints at the underlying production mechanism. However, our understanding of the subject 
is still very poor partially because of the reason that in soft hadronic or nuclear collisions, 
the role of perturbative QCD is less clear. Several new experimental information on mul- 
tiparticle production have accumulated in recent years. Consequently several theoretical 
approaches have been developed in an attempt to understand or organise the data. Sev- 
eral efforts [3-8], which are based on hydrodynamic descriptions, have also been proposed 
to explain the enormous experimental data coming from a variety of collision experiments 
ranging from very low energy (e.g., SIS) to very high energy (e.g., LHC). In addition to 
this, several approaches based on statistical production of particles have recently appeared 
in the literature [9-23]. The statistical models are based on the assumption of local filling 
of available phase space according to statistical laws. Furthermore, there are models which 
try to explain the rather complicated heavy ion colliding processes by analyzing them in 
terms of the basic quark-gluon interaction processes e.g., additive quark model (AQM) [24], 
dual parton model (DPM) [25] or colour neutralization model (CNM) [26] etc. The main 



difference in these models lies in the assumption about the potential number of participant 
partons which can interact independently. Traditional AQM treats the elastic scattering of 
two hadrons at high energy in terms of Pomeron exchange between two quarks involving 
one from each hadron i.e., the number of participants is limited to the valance quarks only 
[24, 27]. However, in DPM and CNM, the number of participant quarks is unlimited be- 
cause these can involve the sea quarks as well. Many attempts have also been made in the 
past for searching certain systematics or scaling relations which are universal to all types 
of reactions, i.e., lepton-hadron, hadron-hadron, hadron- nucleus and nucleus-nucleus colli- 
sions. Any distinct deviation from these relations observed in the ultrarelativistic nuclear 
collisions will be a potential indicator of a new and exotic phenomenon occuring there [1]. 
In order to find out whether deviations arise from the presence of some exotic phenomena, 
e.g., formation of QGP, where we must know the background arising from nucleus- nucleus 
collisions without the presence of any such phase transition [1]. In this paper, we propose 
a model for the multiparticle production and search hints for deviations signaling an exotic 
phenomenon by comparing our results with the experimental data coming from a variety of 
heavy-ion collision experiments. Our model involves a main assumption that the nucleus- 
nucleus collision can be considered as a superposition of independent quark-quark collisions, 
the number of which is determind by the geometry of the collisions. 

Bascially in QCD, the interaction mechanism between target and projectile nuclei can be 
described as follows : A projectile quark exchanges a gluon with a target quark and colour 
forces are thus stretched between them as well as other constituents because they try to 
restore the colour singlet behaviour. When two quarks thus separate, the colour force builds 
up a field between them and as the energy in the colour field increases, the colour tubes 
break up into hadrons and quark-antiquark pairs are created. In this paper, we have basically 
used the model given by Singh et. al. [28] and obtain the charged hadron multiplicities for 
p — p, h — A, and A — A collisions. This model considers a multiple collision scheme in 
which a valance quark of the incident nucleon suffers one or more inelastic collisions with 
a valance quark of the target nucleon. The quark thus loses energy and momenta and 
produces hadrons in each quark-quark collision. However, all the collisions are independent 
and their effects should be incoherently superimposed. Further in A — A collisions, one has 
to adequately incorporate the number of participating quarks and mean number of quark- 
quark collisions. The rest of the paper is organized as follows : In section II, we will give 



the detailed description of the model used in this paper. In section III, we give the results 
and the discussion and finally in section IV, we will give the conclusions based on this work 
and the future prospects of this study. 

II. DESCRIPTION OF MODEL 

If we assume a universal mechanism of charged particle production in the hadron-hadron, 
hadron-nucleus and nucleus-nucleus collisions, it must be driven by the available amount of 
energy involved for the secondary production, and it must depend on the mean number 
of participant quarks. The main ingredients of our model is based on a phenomenological 
model proposed earlier by Singh et. al. [28]. Let us attempt first to understand the basic 
particle production mechanism in p — p collisions. There are several papers emphasizing 
various types of fitting parametrizations in order to provide a unified description for the 
produced charged particles in p — p collisions at various energies. We have taken some of 
these functions with the values of their constants from the literature and fitted them with the 
experimental data. All these parametrizations fit the experimental data for p— p collisions at 
intermediate and higher energies but usually show some disagreement with the experimantal 
data at lower energies ( see Fig. 1). Therefore, we propose here a new parametrization to 
accomodate entire p — p experimental data [29-38] from low energies to very high energies 
(i.e., from 6.15 GeV upto 7 TeV), in a unified way as follows : 



In Eq. (1), a is the leading particle effect and y/s^ is the available center-of-mass energy i.e., 
y/s^ = y/s — rriB — ttit, where m# is the mass of projectile and mr the mass of the target 
nucleon, respectively; a', b' , d and d' are constants. We find the values a' = 1.8, b' = 0.37, 
d = 0.43 and d! = 0.04 as derived from the best fit to the data [29-38]. 

Now, we can extrapolate this expression for the produced charged particles in the hadron- 
nucleus interaction based on two basic assumptions. First assumption is that the number of 
constituent quarks which participated in hadron-nucleus collisions share the total available 
center-of-mass energy and thus the energy available to each interacting quark becomes 
yfsA/N q in h — A collisions. Now, since s a is the available center-of-mass energy in h — p 
collision for one effective collision of quarks of hadronic beam in the target, therefore, second 




(1) 



assumption is that the total available square of center-of-mass energy s A in h — A case 
becomes u q s a provided each quark suffers on average v q collisions. Thus, the expression for 
average charged hadron multiplicity in h — A collisions can be expressed as follows [28] : 



< n ch > hA = N q 



a + b'ln 



+ c'ln 2 



+ d'ln 3 



a. 



(2) 



In the above Eq. (2), y/sj = (ugSa) 1 ^ 2 , where v q is the mean number of inelastic quark 
collisions in the target nucleus and is defined as : v q — Aa q n N /a q n A . Here A is the mass 
number of the target nucleus and o l qN is the quark-nucleon inelastic interaction cross-section 
and equal to one-third of nucleon-nucleon inelastic cross-section (cr™ N ) i.e., a l q n N « (l/2>)a™ N 
[17]. Also a % q n A ,the quark-nucleus inelastic interaction cross-section and is obtained from 
Glauber's approximation by using the following expression [18]: 
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where profile function D A {b) is related to nuclear density, p(b, z) by the relation : 



D A {b) 



p(b, z)dz. 



We use the following expression to calculate the nuclear density : 
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where R and a are constants for any nucleus having mass number A, p is the normalization 
constant and b is the impact parameter. We have taken the values of the constants from 
the Ref. [39]. Furthermore, N q in Eq. (3), is the mean number of participant quarks and is 
defined as follows : 
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where iV c is the number of valance quarks in the nucleus A. 

The generalization of the above picture for the case of nucleus-nucleus collisions goes along 
the same line and can be given as follows : 
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where ^Jsab = A{y q B 's a ) 1 ^ 2 and the mean number of inelastic quark collision v q B can be 



given as follows : 
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Furthermore, mean number of participating quarks N^ B can be calculated by generalizing 
Eq. (6) in the following manner: 

-N B af A N A a% 



N* B = - 



AB U AB 



(9) 



where o™ B is the inelastic cross-section for nucleus A- nucleus B collision and can be ex- 



pressed in the following manner [18]: 



-in _ r 2 
a AB — 7ir 



^1/3 + B l/Z _ 



1 2 



(10) 
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where c is a constant and has a value 4.45 for nucleus-nucleus collisions. One can see that the 
parametrization as given by Eq. (7) gives the most general relation relating nucleus-nucleus 
collisions to hadron-nucleus and hadron-proton collisions and the values of the parameters 
a', b' , c', and d! remain unaltered which shows the similarity in the role of basic quark-gluon 
interaction in all these processess. In a search for creating quark gluon plasma (QGP), 
greater emphasis is laid on the central or head-on collisions of two nuclei. In such a case, 
we assume that all the quarks of the beam nucleus are wounded and the resulting mean 
multiplicity can be obtained by using Eq. (7) in the following manner : 

< n ch >%g ral = A [a' + b'ln(vfs a ) 1 ' 2 + dln 2 {vf s a fl 2 + d'ln\vf B s a f/ 2 - a] . (11) 

The pseudorapidity distribution of charged particles is another very important quantity 
in the studies of particle production mechanism in high energy h — h and A — A collisions. 
It has been pointed out that (dn c h/dr]) can be used to obtain the temperature (T) and 
density (p) of the QGP [40-42]. Moreover, mid-rapidity density reflects the different stages 
of the reaction. To calculate the pseudorapidity density of charged hadrons, we first fit 
the experimental data of (dn c h/ 'dr})^ =Q by using similar parametrization as was used in Eq. 
(1) but with slightly different values of constants. Further extrapolation of two component 
model [43-45] gives (dn ch /drj)^ as follows : 

&)Z = &Yj (l - x)l, " +xN " vn ' (12) 

where x quantifies the relative contribution from hard and soft processes. The fraction 
x corresponds to the hard processes and the remaining fraction (1 — x) arising from soft 
processes. Thus the contribution from hard processes is propotional to the total number of 
collisions i.e., N^ B u^ B and the portion from soft processes is propotional to the participating 
quark constituents i.e., Nf B . 




FIG. 1: Variation of total multiplicities of charged hadrons in p — p collision at various y / s]vjv- 

III. RESULTS AND DISCUSSIONS 

In Fig. (1), we present the inelastic (filled symbols) and non-single diffractive (NSD) 
data (open symbols) of charged hadron multiplicity in full phase space for p — p collisions at 
various center-of-mass energies from different experiments e.g., ISR, UA5 and E735 [29-38]. 
We use inelastic data at very low energies (filled symbols) because there is no NSD data 
available for these energies and also the trend shows that the difference between inelastc 
and NSD data will be very small at lower energies. Further, we fit this data set with four 
different functional forms. The short-dashed line has the functional form as : a + bs 1 ^ 4 which 
is actually inspired by the Fermi-Landau model [46, 47]. It provides a resonable fit to the 
data at higher ^/stvjv with a = 5.774 and b = 0.948 [29]. However, since a summarizes 
the leading particle effect, it should not be much larger than two. The dotted line has the 
functional form as : a + b Ins + c ln 2 s and it fits the data well at higher ^snn but shows 
a disagreement with the experimental data at lower center-of-mass energies. The dashed- 
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FIG. 2: Variation of total mean multiplicity of charged hadrons in central Au — Au collisions at 
various ^/SjVAr- 

dotted line represents the form a + bs n and it also provides a good qualitative description 
of the data with a = 0, b = 3.102 and n = 0.178 [29]. The solid line represents our 
parametrization given by Eq. (1) and it undoubtedly represents the most reasonable good 
fit to the data starting from very low upto very high v /s/vat. 

Fig. 2, shows the variation of mean multiplicity of charged hadrons produced in central 
Au — Au collision with respect to ^stvjv- We also compare our model results (solid line) and 
the modified Glauber model results (dash-dotted line) [48] with the experimental data of 
AGS and RHIC [49-52]. We would like to mention here that the main difference in modified 
Glauber model with respect to the standard Glauber calculations is that for each nucleon- 
nucleon collision modified model uses the value of af^ei a ^ ^ e corresponding centre-of-mass 
nucleon-nucleon collision energy. We find that the results obtained from our model give a 
good description to the experimental data in comparison to the modified Glauber model 
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FIG. 3: Variation of total mean multiplicity of charged hadrons in central collisions for different 
colliding nuclei at various s/snn- 



TABLE I: Total mean multiplicity of charged hadrons produced in d — Au collision at ^Jsnn = 200 
GeV. 





Colliding Nuclei (n ch ) mmMas 


(n ch ) central 


Our Model 


d-Au 82 


162 


PHOBOS [54] 


d-Au 87+1 


W7±\i 



predictions, especially at lower center-of-mass energies. 

In Fig. 3, we show the variation of mean- multiplicity calculated from our model for 
Au — Au, Pb—Pb and Cu — Cu central collisions with respect to ^/snn- We also compare our 
model results with the experimantal data of RHIC and SPS experiments [49-53]. One can see 
from Fig. 3, that the model results give excellent fit to the separate experimental data points 
for Cu — Cu, Au — Au and Pb — Pb, respectively. This concludes that our phenomenological 




FIG. 4: Variation of pseudorapidity density at mid-rapidity for p — p collisionat various -y/J/vJV- 
Filled symbols are data from inelastic p — p events [32, 53-54] and open symbols are experimental 
data for NSD events [33-36, 58-59]. Solid line is the outcome of our parametrization. 

model works well because it considers the picture of incoherent superposition of basic quark- 
quark collisions in describing the nucleus-nucleus collisions. 

In Table I, we have shown the mean multiplicity of charged hadrons in minimum-bias as 
well as in central d — Au collision at y/s NN = 200 GeV. We have also shown the experimental 
results [54] for this collision for comparison. The good agreement between the experimental 
and our results shows that our model can be used for asymmetric colliding nuclei as well. 

In Fig. 4, we present the inelastic (filled symbols) and non-single diffractive (NSD) data 
(open symbols) of (dn c h/ drj)^ for p — p collisions at various center of mass energies from 
different experiments e.g., ISR, UA5, E735, RHIC and LHC [32, 33-36, 55-59]. Further, we 
fit (solid line) this data set with our functional form of Eq. (1) and obtained the values of 
constants as a' = 1.78, b' = 0.025, c' = 0.041 and d! = 0.0017. The central rapidity hadron 
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FIG. 5: Variation of pseudorapidity density at mid-rapidity for A — A collisionat various ^snn. 
Symbols are data from various heavy ion collision experiments [49-53, 60-61]. 



density thus continues to increase with the colliding energy at LHC in a triple logarithmic 
form. 

In Fig. 5, we plot the variation of (dn c h/dr]) v=0 obtained in our model for Au — Au, 
Pb — Pb and Cu — Cu central collisions with respect to ^/snn- We also compare our model 
results with the experimental data of RHIC, SPS and LHC experiments [49-53, 60-61]. Our 
model results compare well with the experimental data from very low energy experiments to 
very high energy experiments. As we didn't take any effect of final-state interactions in our 
model, it verifies that the hadron multiplicity in A — A collisions is mainly driven by initial 
parton production and the effect of final-state interaction is negligibly small as suggested in 
some earlier models [62-63]. 

In Table. II, we give the predictions for the mean multiplicity of charged hadrons, produced 
in minimum-bias events as well as in central events, based on our model calculations at LHC 



TABLE II: Predictions of total mean multiplicity of charged hadrons produced in heavy ion collision 
experiments. 



Experiment Colliding 


; Nuclei 


(GeV) 


1 \min.bias 
V'ch/cur model 


/„ \ central 1 \ central 
V'ch/our model V L ch/ Re f. [66-67] 


/„ \ central 
\'hv/HSD [64-65] 


LHC 


Pb- 


Pb 


2760 


8895 


14377 15000 ± 1000 




CBM 


Au - 


Au 


3.97 (E lab « 8A) 


328.75 


380.8 


308 


55 


Au - 


Au 


8.1 (E lab « 35A) 


670.8 


837 


650 



and CBM energies. Specifically, we calculate the charged hadron multiplicity in Pb — Pb 
collisions at y/s NN = 2.76 TeV for LHC and the mean multiplicity of charged hadrons in 
Au — Au collisions at E iab = 8 AGeV and 35 AGeV for future CBM experiment. We also 
compare our model results with the other model predictions like hadron string dynamics 
(HSD) model predictions for CBM experiment [64-65] and model prediction of W. Busza for 
LHC experiment [66-67]. 

IV. CONCLUSIONS AND FUTURE PROSPECTS 

We have given a parametrization which correctly describes the multiplicity distributions 
in particle and nuclear collisions from a few GeV upto the highest LHC energies. It involves 
the interactions of the constituent quarks of the colliding beams or beam and target. The 
collision of the two Lorentz-contracted particles results in a complete thermalized system. 
We assume that the production of secondary particles is proportional to the fraction of the 
available energy for the participant quarks. This picture describes consistently the 'soft' 
hadron production in p — p, p — p, A — A and p — A collisions. In heavy-ion collisions, more 
than one quark per nucleon interacts due to the large size of the nucleus as well as due to a 
large travel path to be travelled inside the nucleus [68]. For more central collisions, larger 
number of interactions occur and hence we get larger energy available for the secondary 
particle production. In the most central nucleus-nucleus collisions, all three constituent 
quarks from each nucleon will interact simultaneously and they deposit their energy into 
the thermalized collisions volume. We find that our parametrization gives the excellent fit 



to the p — p data for the entire energy range. We also find that a suitable extension of 
this description works very well for the A — A collisions again for the entire energy range. 
We have also compared our model results with the results obtained from modified Glauber 
model to demonstrate the difference between nucleon-nucleon interaction and quark-quark 
interaction pictures. We have shown the charged particle multiplicity in central as well as in 
minimum-bias d — Au collisions at ^s^n = 200 GeV and we find it compares well with the 
experimental result. This exercise shows that our model works well for assymetric colliding 
nuclei case also. We have also proposed an extension of the two-component model based on 
our model and attempted to get the rapidity distribution and we find that it again works 
well. We have given our prediction for LHC and CBM data as well and compare them with 
some recent model predictions [64-67]. 

In conclusion, we have attempted to draw certain universal mechanism for 'soft' 
hadron production in hadron-hadron, hadron-nucleus and nucleus-nucleus collisions. We 
emphasized that QCD does not work for these cases. So lacking a workable theory, 
we have relied on the phenomenology and all the salient features of our model fit very 
well with the experimental data. We hope that the work done here will throw enough 
light on the multiparticle production mechanism in nucleus-nucleus collisions. This 
also reveals that any deviations observed in the data from the predictions will help 
us in identifying QGP formation. Furthemore, we expect that this model can also be 
extended to accomodate the rapidity, pseudorapidity as well as transverse mass distribution 
of charged hadrons in the whole rapidity region which will be studied in a future publication. 
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